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Abstract. We suggest some principal ideas on origin, statistical properties and 
theoretical description of 1 /f- noise exactly as they for the first time were expounded 
in our preprint published in Russian in 1982, and supplement them with short today's 
comments and selected references, with wish to support improvements of present 
generally poor ideologic and mathematical base of the 1 / f- noise theory. 



"Many things are incomprehensible not because our notions are poor 
but because these things are not in a frame of our notions" 
(Koz'ma Prutkov) 

PREFACE 

As far as we know, the present state of the 1/f- noise problem has no principal dif- 
ferences from that thirty years ago when our preprint had appeared (Preprint NIRFI 
No. 157, published in 1982 by the Scientific- Research Radio-Physics Institute in Nijnii 
Novgorod in Russian Federation). Of course, this is said not about varieties of ob- 
jects under experimental investigation but about conceptual level of their theoretical 
interpretation. Reading of today's scientific literature shows us that the interpreters 
as before are entangled in prejudices cultivated by more than centennial history of 
hypotheses, assumptions and approximations in statistical physics and kinetics. We in 
our preprint just made first attempt to disentangle. Therefore we think that it still 
may be useful for interested readers, all the more that it affected all our later works 
(first of all [2211231 HQ [25]). 

Later, we have discovered N. Krylov's book (26] which is devoted to disclosure of the 
prejudices and gave us principal justification of our own findings. In turn, our own 
investigations of "molecular Brownian motion" ([27] - [32], [S]> [IB]) and "electron 
Brownian motion" ([28], [ID] - [13J, [H], [H]) in statistical mechanics gave confirmation 
of Krylov's ideas. As the result, in particular, now we have at our disposal more 
correct and formally substantiated terminology than thirty years ago. Nevertheless, 
our translation of the preprint is as much "one-to-one" as possible. Some necessary 
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corrections and additions can be found in the mentioned references and in comments 
below designed also as references but with mark / . 
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Abstract 

It is shown that thermodynamically equilibrium Brownian motion of charge carriers 
always possesses low-frequency fluctuations in cliff usivities of the carriers and thus in 
power spectrum of "white" noise of the conducting media, with correlation function 
of the fluctuations decaying by logarithmic law (In t/r ) _1 while their spectrum is 
of 1/oj type. Similar fluctuations are peculiar to electric conductivity and current 
in non-equilibrium state (under external field). It is shown that these fluctuations 
are not related to some slow processes (and macroscopic relaxation times). Statistical 
characteristics of the 1/f-noise are completely determined by microscopic parameters 
of "fast" random motion of the carriers. The 1/f-type spectrum reflects absence of 
long-living correlations at this random motion. 

An exhaustive information about spectrum and magnitude of the 1/f-noise is con- 
tained in fourth-order cumulant function of equilibrium current fluctuations. 

The suggested theory estimates the 1/f-noise intensity in agreement with experimen- 
tal data and explains origin of the empiric "Hooge constant" . 



1. Introduction 

1.1. Experimental data and empiric relations. Electric charge transport in vari- 
ous media and systems is accompanied by characteristic low- frequency noise, - so called 
excess noise, or flicker noise, - which is known almost as long as usual white noise, but 
still has no theoretical explanation. 

Surprising peculiarity of flicker noise is that its spectral power density is fast increas- 
ing as frequency decreases, approximately by law u~ a , and does not show tendency to 
saturation down to minimal measurable frequencies ~ 10~ 6 -j- 10~ 7 Hz. In most cases 
the exponent a is close to unit, a ~ 1 , and then it is said about 1/f-noise. 

Intensity of flicker fluctuations of electric current J(t) or voltage in weakly non- 
equilibrium states (in Ohmic regime) is proportional to squared mean current, therefore 
flicker noise is usually considered as result of fluctuations of resistance or conductance. 
Main relationships of electric current 1/f-noise in homogeneous conducting media are 
reflected by approximate empiric "Hooge formula" (for details see reviews [1-2]): 

Here Sj(u) is spectral power density of current fluctuations, J is mean current, N 
is number of charge carriers in a sample of media, and a is dimensionless quantity. 
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Formula fll.ip reflects one more surprising property of 1/f-noise: indifference of shape 
of its spectrum to system's geometry and sizes, the so-called "zero-dimensionality" of 
1/f-noise. 

Formula f ll.ip usually gives satisfactory description of 1/f-noise in semiconductors, 
solid and liquid metals, electrolytes [1]. In case of semiconductors the quantity a is 
almost independent on temperature T and the number iV . The latter observation 
is evidence of statistical independence of contributions from particular charge carriers 
into full 1/f-noise. Interestingly, in various intrinsic (weakly doped) semiconductors a 
has nearly same order of magnitude, a ~ 0.001 (for the first time this fact was noticed 
by Hooge [3]). At high enough temperatures, values a ~ 0.001 -r- 0.01 characterize 
also metals, though there a is temperature-dependent [1]. For electrolytes also one 
can find a ~ 0.001 /[49]. 

There are three sorts of significant differences of 1/f-noise level from that responding 
to a ~ 10~ 3 . In strongly doped semiconductors the noise is much weaker. There, 
according to Vandamme and Hooge, a ~ 10~ 3 (fi/f^o) 2 , where \i is mobility of carriers 
and /io their mobility in pure material. In inhomogeneous systems (among which, 
seemingly, one should rank also very thin films and wires) the noise may be much 
greater. In metals at comparatively low temperatures one can observe "temperature 
1/f-noise" [1,2,4,5]. It is characterized by a s dependence on conductivity's temperature 
coefficient and sample's thermal contact with surroundings. Perhaps, as was noticed 
in [5], in metals the noise represents superposition of such "temperature 1/f-noise" and 
the above mentioned noise which dominates at higher T (T > 150° K). 

In general, the picture of flicker noise in metals looks much more complicated than in 
semiconductors. In many cases measurements of the exponent a show values visually 
different from one ( 0.8 < a < 1.2 ), so that formula (11 .11) appears too rough. 

1.2. Correlation experiments. "Zero-dimensionality" of 1/f noise. "Zero - 
dimensionality" of 1/f-noise expressively manifests itself in so-called correlation exper- 
iments: small neighbor regions of same conducting sample produce uncorrelated con- 
tributions to its 1/f-noise. This seems striking, since we speak about extremely slow, 
in microscopic time scale, fluctuations. Indeed, 1/f- type spectrum can be expanded 
into sum of Lorentzians, 
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were, as was mentioned, 2t[/VLq > 10 6 s. And we may assume that there is some real 
fluctuation mode beyond either Lorentzian, with relaxation time 2ir/Q. However, it 
is hard to imagine flutuations' mechanism what would be space-local but at that lead 
to time scales up to 10 6 s (if not greater) /[5T]. 

Though, the "temperature 1/f-noise" in metals possesses spatial correlations [4]. 
Therefore Voss and Clarke concluded in [4] that it is caused by mere thermodynamic 
fluctuations of temperature "modulating" conductivity. Now, one can state that this 
conclusion was wrong. The theory of temperature fluctuations, consistent with prin- 
ciples of statistical thermodynamics, does not lead to 1/f-type spectra [1,2,4]. At the 
same time, there are no doubts that the "temperature 1/f-noise" is closely related to 
thermal processes. 

It is easy to accept the complexity of situation with metals. In metals, electrons 
are not only charge carriers but also main heat carriers, so that electric and thermal 
processes are mutually entangled. 

1.3. Thermo dynamically equilibrium 1/f-noise and fourth-order cumulant 
of current fluctuations. The accumulated experimental data give evidence that 1/f- 
noise is by its nature thermodynamically equilibrium [1,2]. This principally important 
circumstance means that unknown processes, what are responsible for 1/f-noise, take 
place first of all in equilibrium state, although do not manifest themselves in correlation 
function and spectrum of electric current or e.m.f. 

Nevertheless, in some things these processes must be reflected even in equilibrium 
case. It is not hard to guess that they must lead to flicker fluctuations of intensity 
of equilibrium white noise, S(t) . This is prompted already by the Nyquist formula: 
S = 2Tg , if one interprets 1/f-noise as consequence of fluctuations of conductance, 
g(t) . Voss and Clarke measured fluctuations of white noise power in thin metal film 
and really found that they have 1 / f-type spectrum [4] . This phenomenon also can be 
termed 1/f-noise (equilibrium now just in literal sense). In out-of-equilibrium system 
under external field it transforms into current (and voltage) fluctuations which will be 
termed below current 1/f-noise. 

It should be underlined that the fluctuating power S(t) , in contrast to J(t) , is espe- 
cially phenomenological characteristics of noise. While J(t) always can be written as 
a function of microscopic dynamical variables of a system, S(t) can not be expressed 
by dynamical language neither through these variables nor through fluctuations of 
thermodynamical quantities: temperature, chemical potential, etc., - characterizing 
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system's quasi-equilibrium states. The matter is that the power (power spectral den- 
sity) of white noise represents kinetic quantity. Its definition should involve statistical 
averaging over ensemble and, besides, integration (averaging) over time. But in such 
way one can rigorously introduce only mean valueu 



while fluctuations of S(t) have no strictly definite dynamical sense. The said concerns 
also fluctuations of conductance g(t) and other kinetic quantities. 

Then, how we can rigorously describe fluctuations of the power and conductance? 
It is very simple issue. Since ( S(t) ) is connected to quadratic forms, in respect 
to current or voltage, fluctuations of S(t) are corresponded by current's statistical 
moments of fourth order (and higher orders). 

Hence, it is necessary to consider fourth moment of current's fluctuations: 



For it always there is a rigorous dynamical expression. The brackets with commas 
inside on the right in (jl.2p mean fourth-order cumulant of current. Equality in ( II. 2p is 
known general cumulant expansion of fourth moment at (J(t))o = . 

We have to underline that information about the power fluctuations and equilibrium 
1/f-noise is hidden in the last term of (II. 2p . i.e. fourth cumulant, which characterizes 
non-Gaussianity of current fluctuations. The equilibrium white noise correlation func- 
tion (J(t) J(t2))o is microscopically fast decaying under increase of \t\ — t 2 . Therefore 
low- frequency processes are reflected by only last term of (11.21) . 

From here an important consequence does follow, that in statistical description of 1/f 
noise it is very necessarily to take into account non-Gaussianity of current fluctuations 
(even if it is very weak in one or another case)§. Of course, not current only, but also 
any other physical random process always is more or less non-Gaussian. Here we meet 
the situation when this statement is of principal importance. 

1 Subscript "0" at angle brackets indicates that averaging is made over equilibrium ensemble. 

2 Notice that, for instance, Gaussian noise with randomly varying intensity becomes non- 
Gaussian random process. For more about analysis of noise non-Gaussianity see also: M.Nelkin 
and A.M.S.Tremblay, J. Stat. Phys. 25 253 (1981). 




( J{t) J(0) ) dt , 



( J(ti) J(t 2 ) J(t 3 ) J(U) ) = ( J(fi) J(t 2 ) >o ( J(t 3 ) J(U) ) + 

+ ( J(h) J(t 3 ) >o ( J(t 2 ) J(U) ) + ( J(h) J(U) ) ( J(t 2 ) J(t 3 ) ) + 
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1.4. Fluctuations of mobility of carriers. A number of authors, basing on analysis 
of experiments, have come to conclusion that primary source of 1/f-noise is fluctua- 
tions of mobilities of carriers [1,2,6]. The Hooge formula appears if one assumes that 
particular mobilities fluctuate independently one on another with spectrum 



where /i is mean value of mobility, a 10 -3 . This model allows also empirical descrip- 
tion of 1/f-noise in non-uniform structures: in various contacts, p — n— junctions and 
other semiconductor devices (and, seemingly, in electron emission), moreover, even in 
non-Ohmic regime [1]. There, as one can think, 1/f-noise arises because of fluctuations 
of carriers' flow to a structural change in turn caused by mobility fluctuations. 

In applications to semiconductors, the hypothesis of mobility fluctuations in many 
cases meets competitive hypothesis about fluctuations in number of carriers due to 
slow tunnel transitions to near-surface states and back [1,2,7]. In this model it is easy 
to obtain 1/f-type spectrum as the sum of Lorentzians (thogh saturating at very low 
frequencies), but it is hard to obtain numeric estimates. This model contradicts to 
observed bulk character of 1/f-noise. Besides, special experiments with hot carriers 
also give evidences in favor of hypothesis of mobility fluctuations [6]. 

Then, natural question about physical origin of the mobility fluctuations does arise. 
If ascribing their ground to some slow processes in crystal lattice (e.g. fluctuations 
in number of phonons / 152] . as proposed in [1]), one would have to expect mobility 
fluctuations of different carriers to be correlated one with another. But experiments 
say about the opposite. If, however, we state absence of inter-carrier correlations, 
then the known "zero-dimensionality" paradox arises: why it is possible to observe 
flicker correlations much longer than mean time of residence of charge carriers in a 
small sample? (It seems as if carrier leaves a sample but nevertheless its correlations, 
associated with 1/f noise, stay there.) Up to now in the literature none mechanism 
of the mobility fluctuations was suggested. Evidently, the mentioned paradox strongly 
complicates the problem /[53J. 

1.5. Principal statements of present work. 1/f-noise as result of absence 
of long-living correlations. In spite of many-year experimental investigations and 
many attempts of theoretical explanation, 1/f-noise still remains mysterious [1,2]. It 
is all the more strange in view of that in other respects the systems under question do 
not display mysterious effects what could be associated with 1/f-noise. 





Various investigators' efforts are invariably directed either to search of some new 
"slow" physical mechanisms, which would lead to wide set of very large relaxation 
times (or correlation times, "life-times", etc.) and flicker spectrum of conductance 
fluctuations, or to "approbation" in this sense of already known mechanisms, such as, 
for example, fluctuations of temperature and carriers concentration (see e.g. [2,4,8-11]). 
This way does not lead to success. 

In the present work, exactly opposite approach to the 1/f-noise problem is suggested 
and substantiated, allowing to explain this physical phenomenon. It is shown that 1/f- 
noise can be connected just to absence of macroscopically large relaxation times and 
caused not by specifically slow processes but by Brownian motion of charge carriers 
in itself, i.e. the same "fast" microscopic processes only which produce diffusion and 
white noise. 

We show that Brownian motion (diffusion) of carriers always is accompanied by 
flicker fluctuations of diffusivity and mobility (white noise power and conductance), 
which naturally and inevitably arise in the course of diffusion irrespective to its concrete 
microscopic mechanism, without any slow perturbations (as e.g. random changes of 
thermodynamical conditions of diffusion). The theory to be expounded below not only 
yields 1/f-type spectrum but also ensures correct estimate of level of 1/f-noise and 
explains magnitude of the empirical "Hooge constant" , a ~ 10~ 3 . 

The key, principal, idea of our approach is that 1 / f-type spectrum is not consequence 
of real long-living correlations, but, oppositely, results from absence of such correla- 
tions, indifference of system to random deviations of "rate" of carrier's diffusion from 
its mean regime. Such deviation is not suppressed by backmoving forces, since the only 
its consequence is mere spatial displacement of carrier, that is system's transition into 
a state which is identical in thermodynamical sense to initial state. 

Just in such way, likely, one can explain the observed (see above) mobility fluctua- 
tions. The formulated idea brings also solution of the "zero- dimensionality paradox". 
Since in reality carrier motion is free of any long-living correlations (has no long mem- 
ory about the past), there is no need to speak about breaking of such correlations§ 
under replacements of one carriers in small "noising" sample by others (by similar 

3 It is clear that, very generally, indifference of a system to some spontaneous deviations and their 
accumulation with time (Brownian displacement of a carrier, in our case) can be formally treated 
both as absence of correlations and presence of infinitely long correlations, although physically the 
first of these statements is true. Formal (but not essential!) analogy is given by a superconductor 
with non-decaying currents. 
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reasons, the carriers' generation-recombination processes do not destroy the 1/f- noise) 



2. Brownian motion and equilibrium 1/f-noise. Fluctuations of white 



2.1. Phenomenological description of current fluctuations. Let us consider a 
simple physical situation as follows. Taking a sample of conducting medium, let us 
short-circuit it by closing it into ring. We are interested in thermodynamically equilib- 
rium fluctuations of electric current J(t) in this closed circuit. We shall assume that 
the medium is statistically homogeneous and that charge carriers move statistically 
independently one on another (which usually agrees with real situation in semiconduc- 
tors). Then it is sufficient to consider random walk of one separate carrier along the 
ring-like circuit. 

Introduce designation v(t) for random velocity of carrier in the ring's direction (thus 
v(t) will be scalar), and 



for its path during time t. Notice that according to (12. ip r(t) counts (with plus or 
minus sign) any complete rotation around the circuit. In the same sense we shall treat 
displacement r(t) under charge transfer by not free but localized carriers (hopping 
conductivity) when v(t) is less natural characteristics of the motion. By this r(t) 's 
definition, conveniently, in a stationary state r(t) is random process with uniform in- 
crements (since v(t) is stationary process). In other words, r(t) represents unbounded 
diffusion, although the system under consideration is essentially bounded. In the very 
beginning we would like underline that our results will be quite insensitive to sizes of 
the system (length of the ring) . 

Of course, v(t) always is more or less non-Gaussian random process. The simplest 
usually exploited model of diffusion ignores this obvious physical fact and assumes 
v(t) and r(t) Gaussian. This can be justified when considering carriers concentration 
fluctuations (just they are taken in mind in literature when speaking about "diffusive 
noise", "diffusive mechanisms", etc.) but not in analysis of electric noise caused by 
random wandering of carriers (in our system total number of "noising" carriers does 
not change during diffusion). As already was pointed out, in Gaussian model intensity 
of noise is a priori constant. 



NOISE POWER 




(2.1) 
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Full statistical information about fluctuations of v(t), r(t), J(t) in stationary equi- 
librium state is contained in characteristic functions (CF) 

e t (ik) = (e ife ^) , (2.2) 
e t (m) = (e 4 " Q W) , Q(t)= [ J(t')dt' (2.3) 



JO 

Here k , u are arbitrary probe parameters. It should be underlined that in principle one 
always can concretize rigorous microscopic expressions for these mathematical objects, 
therefore use of CF by itself does not presume any approximations. However, in practice 
one can not do without some assumptions leading to simple enough stochastic model. 
In Gaussian model, as is well known, 

e t (ik) = e' Dk2t , e t (iu) = e-\ SuH , (2.4) 

where D is diffusivity (diffusion coefficient) and S is spectral power density of equi- 
librium current noise at zero frequency. Here it is assumed that t ^> r M , where r /t is 
correlation time of equilibrium noise. More accurate form of ( 12. 4p can be presented 
with the help of functions 

A t (ik) = j ln<e* r <% = \ In Q t (ik) , (2.5) 

Aoo^A;) ee lim \ \n(e ikr ^) , 

t— >oo t 

and similarly for current. Then in Gaussian model 

Aoo^Jfe) = —D k 2 (2.6) 

Probability density distribution of the path r(t) , 

1 f°° 

W t (r) = — e- ikr 6 t (ik)dk , (2.7) 



2tt 



in this model at t ^> r M is 



W t (r) = (AnDt)- 1 / 2 exp ("^) (2-8) 

Now let us consider on-Gaussian model which takes into account fluctuations of in- 
tensity of current noise /[56j. The concept of such fluctuations is meaningful only 
when they are very slow from viewpoint of the scale r M . Therefore their phenomeno- 
logical description requires to use a model treating v(t) like Gaussian white noise with 
random "modulation' of intensity 2D(t) (after whuch v(t) becomes non- Gaussian). 
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In such model by definition 



^expj^ ik{t')v{t')dt'^ ^ = 
exp^- J* D(t')k 2 (t')dt'^ = (2.9) 

exp \-D f k 2 (t') dt 1 + - f [ K D (t' - t") k 2 (t') k 2 (t") dt'dt" + . . 
{.Jo 2 ,/ Jo 



where the first equality corresponds to averaging over white noise at fixed function 
D(t) , the brackets ( . . . ) denote averaging over D{t) 's fluctuations, D = (D(t))' Q , 

K D (t'-t") = (D(t')D(t"))} f — D 2 = (D(t'),D(t")))' Q 

is correlation function of fluctuations of diffusivity D(t) , the dots replace contributions 
to CF from higher-order cumulants of fluctuations D(t) , and k(t) is arbitrary probe 
function. 

On the other hand, there is general exact expansion of logarithm of the CF (12. 9p 
into series over cumulants of the velocityB 

exp j jf ik(t') v(t') rft'j ^ = (2.10) 

I « ^ rt 

eX P 1 / j T / ( U (*0 ' • • • ' V (tn) )o Htl) ■ ■ ■ Htn) dtl... dt r , 

nl Jo 



n=l 



In fact, this is definition of the cumulants [13,14]. In equilibrium {v(t)) = . Equating 
(I2.9p to (12.101) . in view of k(t) 's arbitraryness, one can obtain for fourth-order velocity 
cumulant 

( v(h) , v(t 2 ) , v(t 3 ) , v(U) ) = A5(h - t 2 ) 5(t 3 - U) Koih ~ t 3 ) + (2.11) 
+ 4 5(tt - 1 3 ) 5(t 2 - U) K D {t t - U) + 4 5(h - U) 5(t 2 - t 3 ) K D {t x - t 2 ) 
From here we find (with taking into account that J °° S(t)dt = 1/2 ) 

oo poo 

( v(t) , v(t') , v(t") , v(0) ) dt'dt" = 2K D (t) (2.12) 



'0 JO 

Clearly, in analogous model of fluctuations of current J{t) for its fourth-order cumulant 
appearing in (jl.2p we must obtain 

! ( J(t) , J(f) , J(t") , J(0) ) dt" = -K s (t) , (2.13) 

0^0 z 



Angle bracket with comma-separated factors inside ("Malakhov's cumulant bracket") means joint 
cumulant of these factors. 



12 

where Kg(t) is correlation function of fluctuating power S(t) of equilibrium white 
noise. 

2.2. Non-Gaussian random walk. Characteristic function of displacement. 

Let us consider CF (12.21) . Its logarithm (12.50 for brevity also will be termed CF. 
Expanding CF (12.50 into series over ik we have 

< - x - > / ■ 7 \ 2,717. 

^-^ 2m ! 

m=l 

We took into account that, because of invariance of laws of microscopic motion in re- 
spect to time reversal, equilibrium Brownian motion is spatially symmetric. Therefore 
in (12.141) only even degrees do appear. From the probability theory it is known that 
full number of nonzero terms of the series (12.14p is always infinite. The only exclusion 
is "Gaussian" case when D n (t) = for all n > 3 . 
According to (j2T0 . (l2~5l) . (l2~T0l) . 

D n (t) = \ f ( v(h) , . . . , v(t n ) ) dh . . . dt n (2.15) 



t . .. 

If observation time t ^> r M then D2(t) = 2D, where D is diffusivity. Assume that 
higher correlators (cumulants) of velocity are fast enough decaying at tj\ — > oo . In 
such case in (12. 140 - 02. 150 limits D n = lim^oo D n {t) do exist. The quantities D n give 
values of n -order poly-spectra at zero frequencies, they can be termed also n -order 
"non-Gaussian diffusion coefficients". From (12.51) . (12.140 we have 

OO / ' ] \ 2,771 1 4 

AooCifc) = D 2m = ~Dk 2 + D 4 -- ... (2.16) 

m=l 

If, however, higher-order velocity correlators by some reason are slow decaying under 

separation of time arguments, then the diffusion coefficients D„ ( n > 4 ) may become 

infinite. It means that CF (12.160 is non- analytical function of ik . In such the case it 

is convenient to use integral representation of CF [12]: 

r°° 2D 

Aoo(^) = / (cos kr - 1) — G{r)dr (2.17) 
J-oo r 2 

The multiplier 2D here is extracted by dimensionality reasonings. From (I2.140 - (I2.150 
it follows that G(—r) = G[r) and 



G(r)dr = 1 (2.18) 
In essence, (12.170 is a kind of Fourier transform (taking into account that Aoo(0) = ). 
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In the probability theory the limit transition r^/t — > , analogous to the above con- 
cerned one, is performed in another, more formal, fashion. There t stays finite and 
at that is turned to zero. The result is some "Brownian" process with infinitely 
divisible increments, i.e. with independent increments [12] /[57J. The integral repre- 
sentation like ( I2.17P is termed Levy-Khinchin representation. A fundamental theorem 
of the probability theory states that the kernel G(r) always is non-negative, G(r) > , 
which just means the infinite divisibility. In other respects this function is arbitrary (if 
only the integral in (12.1 7\i is converging). 

The Gaussian diffusion, when D n = for n > 3 , is single peculiar, degenerated, 
case corresponding to kernel G(r) = 5(r) . 

Under "physical" limit transition, when t grows while is fixed, the limit of CF 
(12.1 Op . generally speaking, may not correspond to strictly infinitely divisible distribu- 
tion, and G(r) can be not strictly non-negative. However, in various physical problems 
the increments of r(t) , Q(t) , etc., at large t practically possess asymptotic property 
of infinite divisibility (even in presence of slow decaying, non-integrable, correlations; 
see example in [15]). 

CF (I2.16P contains complete information about large-scale characteristics of the ran- 
dom walk considered in rough macroscopic time scale. Corresponding approximate 
expression for probability distribution of r(t) follows from (12 .7p after replacement 
t (i/c) — > exp (tA QO (i/c)) . At t — > oo , k — )• , a dominating role is played by first term 
of the expansion (I2.16p . that is diffusion is asymptotically Gaussian. 

2.3. Scale invariance of real Brownian motion (r 2 cx t). The Brownian motion 
under consideration represents a physical process realizing in a thermodynamical sys- 
tem and in equilibrium state. Like other thermodynamical phenomena, this process 
must not depend (on scales much greater than characteristic microscopic scales) on de- 
tail structure of microscopic interactions. Consequently, it should possess some spatial- 
temporal scale invariance. A form of this scale invariance is evidently indicated by 
dimensionality of diffusion coefficient D , the macro-parameter determining large-scale 
properties of Brownian motion (and characteristic law of diffusion, {r 2 (t))o = 2Dt). 

Another parameter, what would compete with D in this sense, could appear only due 
to some slow physical processes having significant influence on statistics of diffusion. 

We assume that there are no such slow processes in our system. Then the only 
additional parameters determining (together with D ) complete set of statistical char- 
acteristics of diffusion {D n (t)} are microscopic quantities describing "fast" interactions 
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of carrier with medium. This means that the scale invariance claimed by the "average" 
law of diffusion, (r 2 (t)) = 2Dt , governs the whole statistical picture of Brownian mo- 
tion. In other words, it must look self-similar when spatial scale changes by A times 
while temporal scale by A 2 times (r 2 oc t). 

Mathematically this statement reads, as one cam see from (12.51) . as follows: 

A 2 A A2 , 

at sufficiently large t and small \k\ 

A'Aoo 

Below we shall show that such scale invariance of Brownian motion (at large scales) 
realizes through flicker fluctuations of diffusivity (rate of the motion) with 1/f-type 
spectrum. 



ik 
A 

and 



— = Aoo(ifc) 



(2.19) 



(2.20) 



2.4. 1/f-noise as natural attribute of diffusion. Spontaneous diffusivity fluc- 
tuations. Logarithmically decaying correlations. First, consider relation (12. 20 h . 
Substitution of (" 12 . 1 71) to (I2.20p . after change of variables in the integral we obtain 

AG(Ar) = G(r) 

This functional equation has two solutions: 

G(r) = S(r) , G(r) = A , (2.21) 

\r\ 

where A = const . The first possibility leads to CF (I2.6p . i.e. to ideally Gaussian 
diffusion which has no place in nature. Therefore let us consider the second possibility. 

Choice of the second of expressions ( 12 .2 If) result in divergency of integral in (12. 17ft 
at r — y . This means that the scale invariance can not be perfect, that is it must 
be violated at small (microscopic) scales, which is obvious from physical viewpoint. 
Consequently, we have to cut-off the integrand, for instance, by setting 

G(r) = (2.22) 
|r| + r 

(below it will be seen that details of cut-off procedure are rather insignificant). Thus 
we introduce characteristic spatial micro-scale Tq . It can not be smaller than mean free 
path (in case of quasi- free carriers) or mean step under hopping conductivity /[58J. 
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Inserting (12.221) to (12.171) . we obtain, at k 2 r^ <C 1 (which corresponds to much larger 
scales than r ) , 



A^ik) ~ — Dk 2 A In (2.23) 

r^k z 



/ [59J . Non-analyticity of this function manifests presence of non-integrable long-living 
higher-order correlations of v (t) . At that, however, it is unpleasant that function f |2.22j) 
does not satisfy the necessary condition (12.181) . and therefore expression (12. 23ft does not 
identify quadratic term oc (ik) 2 . This means that the invariance must be destroyed at 
large scales too, as far as observations of the diffusion process take a finite time. 

Hence, we have to go back to CF (12. 5p for finite time and analyse it with the help 
of relation (12.191) . With this purpose, let us use analogue of the representation (12.221) 
as follows: 

A t (ik) = I (cos kr - 1) — G{r,t)dr , (2.24) 



G(r,t)dr = 1 (2.25) 

3 

At t — > oo function G(r, t) should turn into f !2.22|) (but now with coefficient A 
depending on t because of condition (I2.25p . The condition ( I2.25P says that at k — > 
(and t ^> r M ) CF (I2.25P tends to expression —Dk 2 , i.e. is asymptotically Gaussian. 
One can verify that these reasonings together with (I2.19p and (I2.24p imply the following 
form of the kernel in representation (I2.24[) : 

G(r,t) = t^-f(^-\ , (2.26) 



where D' is a constant with same dimensionality as D has, A(t) is determined by 
the normalization condition (I2.25p . function F(z) satisfies requirements 

POO 

F(0) = 1 , F(z) -^0 at z^oo, / F(z) dz = 1 (2.27) 

Jo 

(they always can be satisfied due to presence of free parameters A{t) , D'). With use 
of (I2.27P we find from (I2.25P that, regardless of concrete form of F(z) , 

^ = lf ^^rv.niy 1 (, 28) 



AD't J \r\ + r J \ r o 

at t^> tq . Here a microscopic time scale has appeared: 



T = — (2.29) 
u 2D' K J 
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Next, insert fl2T26|) . ( I2T28]) to (12T241) and consider first terms of expansion (jZ2j]) : 



At(iJfc) = -.DA; 2 + -DD'tA{t) - ... , 
3 

£) 4 (t) = 8 DD'tA(t) (2.30) 

The second term of series (12.301) contains information about diffusivity fluctuations. 
Indeed, phenomenological relations (12. lip . (12. 12[) and definition (I2.15P imply 



d 
dt 



— W 4 (t) = 2AK D {t) , (2.31) 



and then comparison of ( 12.301) and ( 12.311) yields (for t ^> r 

1 „„, d 2 o / t V 1 2 „„, / t 



KMt) -i DD 'lf- t2 {^7 a ) ~l DD '{ la 7 a ) < 2 ' 32 > 

Thus, correlation function of diffusivity fluctuations decays by logarithmic law. It is 
obvious, already from dimensionality considerations, that corresponding spectrum is 
of 1/f-type. 

We see that the flicker fluctuations appear as generic property of real Brownian 
motion (which inevitably becomes scale-invariant when all microscopic scales "are left 
behind") /[£!]. 

Interestingly, these fluctuations rather weakly tell on shape of the displacement's 
probability distribution ( 12. 7ft . As a not complicated analysis of (12.241) . (12. 26ft does 
show, at t ^> r and k 2 r\ <^ 1 CF (I2.24[) has approximately universal form 

AtCifc) = j^ ln (^ 2 + ^) , (2.33) 

TO 

where c is a quantity of order of unit. CF (I2.33f) rather weakly differs from ide- 
ally Gaussian one, (12. 6ft . although coefficients of expansion of (I2.33H into series ( 12.141) 
unboundedly grow with time. Correspondingly, difference between distribution ( 12. 7p . 
resulting from (I2.33p . and Gaussian "bell" ( 12. 8 \ is almost unnoticeable. In other words, 
flicker fluctuations D(t) do not destroy usual picture of diffusion. If we "let out" an 
ensemble of Brownian particles to walk from some point, their distribution law (evo- 
lution of their concentration) is almost Gaussian bell (12. 8p / [62] . Moreover, it can be 
shown that fluctuations D(t) practically have no influence on (equilibrium) concen- 
tration fluctuations, so that the latter can be well described in the frameworks of usual 
"ideally Gaussian" model of diffusion. 
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2.5. Spectrum of diffusivity fluctuations. Quantitative estimates. Origin of 
the "Hooge constant" . Consider spectral power density of (relative) fluctuations of 
diffusivity, 

S D (u) = 2_ 
D 2 ~ L> 2 



Ssd(uj) = -^r 1 = ™ / K D {t) cos utdt (2.34) 



We omit calculation of the Fourier integral. The result, in the region of our interest 
(wr «C 1 ), with good accuracy is equal to 

3D \u)\ ( m \uj\tq ) 2 

Let us compare this expression with the empirical formula (jl.3p . Slow flicker fluctu- 
ations D(t) and fluctuations of mobility, //(£) , should be connected, as it follows from 
simple physical reasonings, by the Einstein relation D(t) = Tfi(t) (in Sec. 3 we shall 
rigorously prove it for steady non-equilibrium state). Consequently, S$d(u)) = Ssfj,(u>) , 
and we can rewrite (12.351) as 



„ , , 2n a(u) 

SsM = — ^ , (2.36) 



\U)\ 



D' 



Spectrum of relative fluctuations of power, S(t) , of summary white noise produced by 
iV statistically independent carriers results from (12.351) after division by iV : 

S ss (u) = Ssd(u)-^ (2.37) 

The parameters r and r are micro-scales representing lower bounds of the scale 
invariance of Brownian motion. Physically, it seems obvious that in homogeneous 
medium 2Dt > , i.e. D' < D . Indeed, if already after single typical "free path" 
(or "hop"), with length Ao , self-correlation of direction of motion vanishes, then r 
must turn to Ao , and therefore quantity 2Dt$ can not be essentially smaller than 
A^r 2 sm- 

If temporal invariance takes shape starting just from minimal accessible scale ~ 
Tq/2D , then 2Dtq w Tq and D' ~ D . In this simplest case tq w r M (with denoting 
typical free path time, or time between hops from one localized state to another). At 
that, "width" of function F(r 2 /Wt) in (12T26|) coincides, in view of d2~27j) . with width 
of Gaussian bell ( 12. 8p . In this remarkable characteristic case, diffusion is invariant to 
maximal extent, since it is described by only two parameters, D and tq . 



2 
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Confining ourselves by this situation|j, let us set 2Dt = r 2 , (D' = D) and estimate 
the "Hooge constant" a(cj) (compare (12.361) with (II. 3p ). 

In semiconductors, typical free path time ~ 10~ 12 s. Taking r ~ 10~ 12 s, at 
frequency u/2tt = 1 Hz we obtain a = (l/3)(ln uutq)~ 2 f» 5-10~ 4 , and a ~ 1.2-10 -3 at 
frequency u/2tt = 10 4 Hz , in good agrrement with typical experimental values ~ 10~ 3 
(see Introduction). Thus our theory gives correct estimate of the Hooge constant and, 
hence, 1/f-noise level in many real situations [1]. 

We see that the "Hooge constant", a ~ (l/3)(ln ur )~ 2 , is determined by micro- 
scopically small time r and total duration ~ 2h/uj of 1/f-noise observations only. 
Frequancy dependence of a(u) is very weak. But it ensures integrability of spectrum 
(I2.35l) -( l2~36l) at low frequency (as it should be for stationary fluctuations) /[6l]. 

If 2Dt > r 2 then, according to (T2T361) . 1/f- noise level essentially decreases. This 
statement agrees with observed lowering of noise in strongly doped semiconductors 
[1]. Indeed, there ro is determined by carriers' scattering by impurities. However, the 
scale ro must be determined by much more slow interaction with phonons. Therefore 
2Dt 3> Tq . Let us assume that 2Dt rs Aq with A being mean free path under 
interaction with phonons only (i.e. in pure material). Then in (12.361) Tq/2Dtq Ri Tq/Aq , 
which can be written also as t 2 /2Dtq rj (/i//io) 2 , in agreement with the empirical 
Vandamme-Hooge formula (see Introduction) / [65] . 




2.6. Microscopic origin of flicker fluctuations. We demonstrated that if Brownian 
motion has no inner macroscopic scales and therefore possesses natural scale invariance, 
r 2 oc t , then it is accompanied by flicker fluctuations of its intensity ( "rate" of diffusion). 
Already the very premise implies the statement from Introduction that related long- 
living correlations are manifestation of absence of long-living memory in random walks 
of carriers, i.e. these correlations are imaginary (nevertheless, one has to describe them 
statistically like real correlations) ^|66j. Both their behavior reflected in ( 12.321) and 
logarithmic cut-off of spectrum (I2.35I) - (I2.37I) at very low frequencies are determined by 
the microscale tq only. In essence, the words about "imaginarity" of flicker fluctuations 
and absence of macroscopic scales are equivalent. 

The "microscopic" interpretation of obtained results enforces to refine some conven- 
tional notions. The usual view at fluctuations of mobility and conductance is that 
these kinetic characteristics always, at any time moment, have exactly definite instant 

5 Expectedly, it takes place, in particular, under domination of one type of inelastic interaction of 
carriers with thermostat, when all correlations of v(t) are fast decaying during time - rl/2D . 
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values which, however, undergo some slow random perturbations. Then one should find 
mechanisms of the perturbations, though such attempts constantly are unsuccessful. 
But there exists quite different point of view: kinetic characteristics have no certain 
"instant" values, and just this physical fact explains 1/f-noise. 

Really, it is not hard to understand that interaction of Brownian particle (carrier) 
with thermostat not only impels it to diffusive motion but simultaneously causes de- 
viations of its random motion from some average regime. The latter can be definitely 
and quickly measured if observing a large ensemble of carriers, as is usually done in 
experiments. However, a separate carrier "knows nothing" about properties of the en- 
semble and is not obligatory at all to display walk with certain diffusivity and mobility. 
The latter concepts have a meaning only in respect to a large enough set of realiza- 
tions of random motion, while their application to individual random (but dynamical) 
trajectory of carrier is senseless. A single separate carrier merely has no instant diffu- 
sivity and mobility, so that one can speak about fluctuations of these kinetic quantities 
in suitable phenomenological language only. Both D , \x and correlation functions 
K D {t) , K^{t) , etc., are characteristics of ensemble of trajectories. 

Of course, if a thing does not exist (as clear dynamical characteristics of motion) 
then it can not fluctuate. Therefore we have to complicate the picture. The diffusion's 
property what is empirically perceived as flicker fluctuations D(t) (or fi(t) ) is mere 
consequence of temporal uniformity of Brownian motion, thermodynamical indifference 
of the system to "rate" of diffusion. Wherever a carrier has occurred at given time mo- 
ment, it each time "starts from beginning", and its past is of no importance. Therefore 
any deviations from ensemble-average regime of motion do not induce a compensation 
and accumulate in time|j • These deviations, or fluctuations of "degree of randomness" 
of motion, cause neither dynamical nor thermodynamical reaction of the system /|71j. 
At the same time, they stay in frames of characteristic diffusive law r 2 oc t . The result, 
as we have seen, is 1/f-noise. 

All that picture, which is slipping away from conventional kinetic's ideology, in prin- 
ciple could be described in formally rigorous and complete enough way in terms of 
the fourth-order correlators (cumulants) ( II. 2p . (" 12 . 12j) . (12 . 13j) . if we were able to ana- 
lyze them by methods of statistical mechanics. However, such the task is extremely 
complicated, even in case of quadratic correlators. 



6 Notice that in [8] an attempt was made to introduce to a model of diffusion similar non-decaying 
("residual") correlations, but they were addressed to fluctuation in concentration of diffusing quantity. 
In our model, long correlations automatically appear in diffusivity fluctuations D(t) . 
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In essence, all the aforesaid concern any kinetic quantities, not "electric" ones only. 
Thus, likely, we can say that in any system, where time-uniform transport of some 
extensive physical quantity takes place, there are flicker fluctuations of spectral power 
of random flows (Langevin forces) and, in non-equilibrium state, also similar fluctu- 
ations of irreversible flows. For example, the heat transfer must be accompanied by 
flicker fluctuations of thermal conductivity and, under temperature gradient, heat flow. 
Another example is given by observed fluctuations (with 1/f spectrum) of energy losses 
in quartz resonators (see e.g. [9]) /[72"]. 

2.7. General description of equilibrium 1/f-noise. Now, let us take arbitrary 
two-terminal dissipative electric conductor whose leads are shortly connected. Consider 
equilibrium diffusion of charge Q(t) (see (12.31) ) around this closed circuit. Here, also 
the usual diffusion law is satisfied on average, {Q 2 (t))o = St . If transport of charge 
is determined by microscopically scaled processes only, the "charge Brownian motion" 
Q(t) again must be scale-invariant. Obviously, now we can immediately write an 
expression for spectrum, Sss(w) , of relative fluctuations of the white noise power. It 
is sufficient to make in (I2.35I) - (I2.36I) replacements 2D =>- S , r =>■ g , where g is 
characteristic microscopic scale of charge transfer, while r has the same meaning as 
before. Thus, at (ur ) 2 <C 1 =,, we find 



If the conductor represents a homogeneous sample with length L , then it is easy 
to transform (I2.38j) back to f!2.35[) - (I2.36j) by setting g = er /L, where e is electron 
charge, and expressing S via diffusion coefficient and number of (independent) carriers. 
Here, go is charge transported through leads in outer circuit (or short-circuit) when 
an inner carrier moves to distance r$ . 

Next, consider, as an example, p — n— junction In this case S = e 2 n , where n 
is mean number of carriers crossing the junction per unit time. If carriers are not 
correlated one with another, then, evidently, we can write go ~ e • Correspondingly, 
at r ~ 10~ 7 s and u/2tt = 1 Hz we obtain from (12.381) estimate 



The scale To must be determined by "life-time" of carrier on the junction, so that 
T n ^> 1 . Relation of r to physical characteristic times of the system can be revealed 
in a more detail model only. 




(2.38) 




(2.39) 
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Of course, in the framework of the presented phenomenological statistical theory we 
can not generally expect estimates better than by an order of magnitude^. 

3. Current 1/f-noise in steady non-equilibrium state and nonlinear 
fluctuation-dissipation relations (fdr) 

3.1. Characteristic function of transported charge. Cubic FDR. The use of 

mathematical formalism of non-Gaussian random processes gave us possibility to con- 
sider 1/f-noise in spite of leaving it in equilibrium state. Under switching on an electric 
field, fluctuations of diffusivities of carriers, D(t) , transform into fluctuations of their 
mobilities, , and thus of conductance, g(t) , and become source of current 1/f- 
noise. For correct analysis of this non-equilibrium state we shall apply the nonlinear 
fluctuation-dissipation relations (FDR) described in [16]-[21] /[73]. 

Let us return to above mentioned tw-terminal conductor. Let before time moment 
t — it is in equilibrium with its surroundings (and short-circuited) but after t = is 
subject to a given constant voltage x(t) = const . In such the shown in [16,17] 

(see also [18,21]), due to the time reversibility of microscopic dynamics, the following 
exact generating FDR takes place [ 



exp 



iu(0 - |] J(t')^j) 



expl- iu{t-t')J{t')dt'\) , (3.1) 



where u(t) is arbitrary probe function. The subscript "x" reminds of non-equilibrium 
character of fluctuations. Taking here u(t) = u = const , introduce CF 

1 f l 

A t (iu\x) = -ln(exp{ra / J(t')dt'}) x (3.2) 
t Jo 

In equilibrium case it coincides with (12. 51) . From ( 13. ip it follows that 

A t (iu — x/T\x) = A t (—iu\x) (3.3) 

Expansion of ( 13 .3[) into series over iu implies an infinite chain of FDR connecting 
average value, {Q{t)) x , of charge transported during observation time t, and various 



7 The authors of [1] quite successfully describe 1/f-noise in non- uniform semiconductor structures 
on the base of empirical model of mobility fluctuations (see Introduction) with a w 0.001 . Interesting 
exception is presented by diodes with the Shottky barrier [1]. 

8 This formula concerns the conductor itself and all the surroundings in contact with it, assuming 
their mutual thermodynamical equilibrium at common temperature T before an external perturbation 
(electric voltage). It is asumed also that there is no external magnetic field. 
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cumulants of fluctuations J{t) and Q(t) . With the help of these FDR, it can be 
shown, in particular, that 

oo „ , 

/ 7» \ Im— 1 

<Q(0>. = £*»(?) <<3 (2m) (*)>*, (3-4) 

m=l 

where 

(Q (n) W),= / (J(t 1 ),...,J(t n )) x dt 1 ...dt n (3.5) 

are n -order cumulants of transported charge, and the numbers defined by 

formulae 

z 1 1 1 

2 ' Cl ~ 2 ' ° 2 ~ ~24 ' ° 3 ~ 240 ' " ' 

m=l 

Differentiation of H3.4[) in respect to t leads to fluctuational representation of average 
current: 

(J(t)) x = ^j\j(t),J(t')) x dt'- (3.6) 

1 /£\ 3 

~ 6 VT 

Consider weakly non-equilibrium state and expand quantities in first row of (13. 6j) 
into series over x : 



^c m z 2m 1 = tanh - . r, = - . r 2 = -— . <v, 



( J(t), J(ti), J(t 2 ), J{t 3 ) ) x dhdhdh + 



(■/(*)>* = gi{t)x + g 3 (t)x 3 + ... , (3.7) 

( J(t), J(t') } x = ( J(t), J(t') ) + x 2 R(t, t') + ... (3.8) 

(thus, for simplicity, we assumed the conductor electrically symmetric). Extracting in 
(13. 6 p linear terms, in respect to x , we find 

<fc(t) = I l\j{t),J{t')) Q dt' , (3.9) 



that is usual fluctuation-dissipation theorem (FDT). Extraction of cubic terms yields 
"cubic FDT": 



1 '* 



<&(*) = y y o R(t,t')df - (3.io) 
( . v .. ( , x j{t),j{t 1 ),j{t 2 ),j(t 3 )) dt 1 dt 2 dt 3 

Notice that similar quadratic and cubic FDR were in part investigated by Efremov [9] 
and Stratonovich [20]. 
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3.2. Weakly non-equilibrium state. Statistical expression for correlation 
function of conductance fluctuations. Let us consider relation (I3.10p at t ^> , 
where r M is correlation time of equilibrium current noise in short-circuited conductor. 
The function gs(t) , representing cubic part of average current response, covers con- 
tributions from both "fast" electric processes (responsible for conductivity in itself) 
and comparatively slow thermal processes, including change of conductance because of 
Joule heat. Hence, generally g 3 (t) , in contrast to gi(t) , is not constant at t ^> but 
can contain slow long "tail". Moreover, this tail may be non-stationary. 

The function R(t,t') , as one can see from (|3.8p . describes "excess" fluctuations of 
current in weakly non-equilibrium state, i.e. in quadratic regime in respect to voltage 
(and mean current). Hence, in particular, R(t,t') contains information about flicker 
fluctuations of current. From (13.101) and ( 12.181) it follows that low-frequency (flicker) 
current fluctuations in the quadratic regime are determined, first, by equilibrium fluc- 
tuations of white noise power and, second, by slow or non-stationary thermal processes 
(caused by Joule heating). Notice that the second source, possibly, is related to the so- 
called "temperature" 1/f-noise what is sometimes observed in metals (see Introduction) 

Here, we are interested in the first source of current 1/f-noise having thermodynami- 
cally equilibrium nature. Therefore we shall neglect non-equilibrium thermal processes, 
i.e. suppose that at t > r M the system tends to strictly stationary state. In such the 
case g%{t) turns to constant at t > r M . Besides, at t ,t' 3> t m the function R(t,t') 
becomes a function of the time difference only: R(t, t') = R(t — if) . Then we can ex- 
ploit he phenomenological concept of conductance fluctuations and write (at t , t' ^> 
and, clearly, \t — t'\ 3> ) 

R(t,t') = R(t-t') = K g (t-t') , (3.11) 

where K g {i — if) is correlation function of these fluctuations§. Next, differentiate 
( 13.1 Op in respect to t . Since at that the left side of (I3.10P vanishes /[76J, from A3. 10|) - 
( l3~TTj) and ( 12T3|) we obtain 



K 9® = zflj I ( J ( t )' J ( t ') ) JH"), J(0) >o dt'dt" , (3.12) 

K 9 (t) = K s (t)/AT 2 (3.13) 



9 Notice that use of model of conductance fluctuations automatically presumes stationarity. In other 
words, contribution of non-stationary processes (if any) to current 1/f-noise can not be described in 
this model surely exploited in the literature. v^|75j 
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Thus, correlation function of conductance fluctuations (caused by equilibrium pro- 
cesses) is expressed through fourth-order cumulant of equilibrium current fluctuations. 
In essence, formula f 1 3 . 1 2 j) must be treated as statistical definition of K g (t) , since the 
fourth equilibrium cumulant of current always can be represented in rigorous terms of 
statistical mechanics. 

Because of (13.121) and the Nyquist formula, S = 2Tg , spectra of relative flicker 
fluctuations g(t) and S(t) are identical and (in weakly non-equilibrium steady state) 
both coincide with spectrum of relative fluctuations of current. From (12.381) and (I3.12p . 
in the quadratic regime, we obtain 

Sj{w) = J" ■ tt^ " o ( ^ Mt )" 2 , (3.14) 

where J = gx . In case of homogeneous conductor in the approximation of independent 
carriers formulae (12.351) . (I2.37P and (13.121) yield 

Sj ^> = J at I I ' a w = 7TT~ ( w r ° 3 - 15 

As was noticed above, this result is in rather satisfactory quantitative agreement 
with general empirical formula ( II. ip . Since we based om "first principles" only, such 
result can be qualified as important evidence in favor of our approach. 

3.3. Nonlinear conductor. Current 1/f-noise in non-Ohmic regime. Our anal- 
ysis is valid for any (nonlinear) conductor, but in the Ohmic (quadratic in respect to 
current) regime. For description of non-Ohmic regime it is necessary to take into ac- 
count non-Gaussianity of white noise in itself which, according to the FDR [17,21], is in 
close relationships with dissipative non-linearity. It is natural to suppose that character 
of these relationships is not affected by slow flicker fluctuations. In the corresponding 
statistical model all kinetic parameters of white noise fluctuate in coordination with 
each other. This can be mathematically formulated (for steady state at t ^> r M ) as 
follows (compare with (I2.9P ): 

1 /"* 

A t {iu\x) = - In (exp { / £(t ') S{iu\x) dt '})' x , (3.16) 
t Jo 

where Jl^{t')dt' is already considered process of scale-invariant diffusivity fluctu- 
ations, with £(t) = S(t)/S , and S{iu\x) is CF of white noise /[77j. The latter, 
according to (I3.3p . should satisfy similar relation 

S(iu — x/T\x) = S(— iu\x) (3-17) 
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Of course, in general S , as well as r and r , is dependent on x . 
Let us write the white noise CF as series 

- (iu) 2 

S(iu\x) = iuJ{x) + ^-S{x) + ... , (3.18) 

where J(x) is average current, and S = S(x) spectral power of noise. Combining 
(13.161) with (13.181) . it is easy to obtain, for current flicker fluctuations, 

KM = l\x)^- fin (3.19) 

and, correspondingly, 

Sj(w) = " 7 ■ -r-r (In \u\To ) (3.20) 

6Tob{X) \0J\ 

Thus, spectrum of current fluctuations is presented, as before, by expression like 
(12.381) . but now with taking into account possible dependences on x . At that, FDR 
(13.171) dictates definite, may be rigid, relations between J(x) and S(x) . 

For example, consider again an ideal semiconductor diode (p — n— junction) in non- 
Ohmic regime, or shot noise regime. In this known, 

S( x ) = e\J(x)\ (3.21) 

at e\x\ ^> T . Assuming that qo = e and tq do not change, we find from (l3.20[) -( |3.2in 
that 

Sj{ U ) ~ ■ (3.22) 

Hence, in nonlinear regime intensity of current flicker noise grows proportionally to 
mean current, that is merely to number of carriers crossing the junction per unit time 
(see [1]) /[7S]. 

4. Conclusion 
Let us resume main aspects and results of the present work. 

1 . The experimental situation around 1/f- noise does not fit in frames of traditional 
notions. The matter is in confusion between concepts concerning statistical properties 
of ensembles and concepts concerning individual charge carrier random walks. The 
confusion is unconsciously provoked by use of traditional Gaussian model of diffusion. 

2. Our principal result is that 1/f-noise can be explained without attraction of 
special physical mechanisms, merely as generic property of random Brownian motion 
of charge carriers, i.e. the well known process involved in most of electric phenomena. 
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3 . We constructed a simple novel model of Brownian motion, exploiting only the 
same physical premises as the standard model do, but refusing assumption about 
Gaussian character of statistics of real Brownian motion, since this is physically sense- 
less idealization /[HQ]. Thus our theory is more adequate (though, of course, may 
require further improvements). 

4 . The constructed theory inevitably leads to conclusion that real Brownian motion 
always possesses fluctuations of "rate of diffusion" with 1/f-type spectrum. 

In the ensemble language one can say that diffusion coefficient and spectral power of 
"white" electric noise undergo fluctuations with 1 / f-type spectrum. In non-equilibrium 
states these fluctuations manifest themselves in mobilities of carers, conductance and 
current. 

5 . Random motion of separate carrier (and, generally, any concrete Brownian trajec- 
tory) has no certain diffusivity (diffusion coefficient). What is for the "rate of diffusion" , 
its spontaneous (thermodynamic) fluctuations have no characteristic time scale, since 
make not back impact upon dynamical picture of motion or thermodynamical state of 
the system. 

1/f-noise is just result of absence of long-living correlations in mechanisms of random 
motion of carriers (and has no relation to some macroscopically large relaxation times). 
Therefore, spectral contents of 1/f-noise is indifferent to system's sizes. 

6 . Intensity of 1/f-noise is determined by microscopic scales only. In the constructed 
model, there are only two such parameters which indicate spatial and temporal lower 
bounds of scale invariance of Brownian motion. 

7. The resulting estimate of 1/f-noise intensity is in satisfactory agreement with 
experiments. At that, origin of the small empirical "Hooge constant" is revealed. 

8 . A relevant microscopic information about 1 / f-noise in stationary states can be 
obtained, in principle, from fourth-order equilibrium correlators (cumulants), with use 
of rigorous methods of statistical mechanics. This is conceptually quite new problem 
of theoretical physics. 

Appendix 

Let us show now in the spectral representation, i.e. in the frequency domain, that 
fluctuations of power of non-Gaussian equilibrium white noise and corresponding con- 
ductance fluctuations possess 1/f-type spectrum. 
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Let J(oS) be Fourier transform of current fluctuations. According to the Wiener- 
Khinchin theorem, 

(J(w),J(a/))o = 2tt S 5 (u + u') , 

where S = const for white noise. From here it follows that 

( J(X 2 lo), J(AV)) = 2ttS6(X 2 (uj + uj')) = 
n S I J (a) J(w') 

Hence, 

A J(X 2 co) oc J(w) , (.1) 

where symbol oc means identity in statistical sense. If current fluctuations have no 
coupling with some slow processes, then the scale invariance expressed by formula Oj) 
must extend to the whole statistics of the fluctuations. Taking this in mind, consider 
fourth-order spectral cumulant: 

( J(wi), J(u 2 ),J(u) 3 ),J(u)±)) = 2irS4(ui,u 2 ,u 3 )6(u) 1 + u 2 + u; 3 + U4) , (.2) 

where S4 is tri-spectrum. According to the analogue of (12.111) for current, 

S 4 (u;i, ^2,^3) = + w 2 ) + S s (ui + u 3 ) + S s (oj 2 + uj 3 ) (.3) 

Replacing in ([72]) J{oj) by A J{\ 2 ui) and combining the result with identity ^TTJ) and 
expression ([72]) itself, it is easy to deduce that 

A 2 ^(AVl, \ 2 u 2 , \ 2 u 3 ) = S^Ux, u 2 , u 3 ) , 
X 2 S s (\ 2 co) = S s (co) 

From here we find 

const 

SsM = —— , (.4) 

that is the power and conductance fluctuations have 1/f-type low- frequency spectrum 
(solution Ss(u>) oc S(u>) is not appropriate since ^4 must be a smooth function). Our 
above results differ from this more formal one only by slight distortion of ideal scale 
invariance. 
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More adequate logics should take into account that slowly decaying statistical correlations in 
general do not presume slowly decaying cause-and-consequence correlations (long- living memory). 
In opposite, long statistical correlations (or even not decaying ones) can arise from indifference 
of system to pre- history of transport processes in it (i.e. just from absence of long memory)! For 
explanations see e.g. [28l [37] and introductory and/or discussion-conclusion sections in other our 
works, [27] - [?]. 

In brief, if charge conduction is indifferent to amount of already transortted charge Q{t) , then it 
is unable to ensure well certain value of time- averaged current Q(t)/t , which just means presence 
of arbitrary long statistical correlations between time-local current values. 

That is what must be expected from rigorous statistical mechanics [35], in place of its crucial 
roughenings in conventional kinetics (pathological simplifications "mutilating" the very way of 
thinking of theoreticians). 

[54] M.B. Weissman, Rev. Mod. Phys., 60, 537 (1988). 

[55] One can not break long correlations in a physical process which by itself is broken into many 
"short" random events! Long statistical correlations between constituent events of such process 
say only that all the events are equally responsible for deviations of the whole proces from some 
"average" regime. The deviations, in turn, result from independence of presently occuring events 
on total number and time distribution of previously happaned ons. 

Such "purely random" process has no certain a priori "probabilities" for constituent events (free 
flights and collisions of carriers, or their passes through a sample, etc.). Therefore an adequate 
theory should start from dynamical (Hamiltonian) model or, at least, deal with full random 
realizations of the process as the whole, while use of a stochastic model with a priori postulated 
partial probabilities automatically "kills" 1/f-noise in theory and makes it a mystery in practice. 

[56] Quantities like D and S , i.e. spectral densities, by their nature are time-nonlocal: they are 
attributes of not definite time moments but time intervals essentially greater than characteristic 
memory life-time of system under consideration. Therefore, they are out of scope (control) of 
the system. Consequently, their fluctuations do not cause system's back reaction and hence are 
scalclcss (indifferent to time averaging). 

Similarly, if you constantly forget your past, then you can not plan your future, and uncertainty 
of results of your activity will be indifferent to its duration. 

Generally, any real (physical) memoryless random process possesses long-living scaleless 
("flicker") fluctuations in its kinetic properties and thus is non-ergodic (terrible word!) in 
respect to them. 

[57] Here, the statistical independence is taken in mind, i.e. independence in the sense of probability 
theory where it by its definition means that joint probability distribution of different incre- 
ments is factored into product of their particular (marginal) distributions. In physics, scientists 
naturally take in mind "independence" in physical cause-and-consequence sense. The great mis- 
take in practical statistical physics of noise, fluctuations and transport phenomena is common 
opinion that the "physical independence" can be identified with statistical one, thus giving rights 
to replace "honest" analysis of many-particle dynamics by 'hand-made" probability-theoretical 
models composed of a priori prescribed probabilities and glued by statistical independency. 
However, Krylov showed that this is wrong opinion, and rigorous statistical mechanics does not 
give such rights |26j . There phase trajectories of many-particle system can not be divided into 
statistically independent "elementary" events. Or, in other words, real events have no definite a 
priori probabilities, while their a posteriori probabilities are different at different phase trajec- 
tories (experiments). 
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In the commented preprint, the r(t) 's increments will loss their statistical independence and 
become all mutually correlated immediately when the kernel G(r) from (|2.17[) will acwuire time 
dependence in (|2.24l) . 

To be more precise, in general there are more than one characteristic "free path" lengths, and 
apparently r can not be smaller than minimum of them. For example (touched in Introduction 
and below), if energy relaxation of carriers is more slow than velocity relaxation, then we have 
at least two different free paths, "elastic" and "inelastic". Then it seems natural to equate r$ 
to shortest of them ("elastic"), while associating \J2Dtq with the longest ("inelastic"), so that 
t'q/2Dtq = D'/D <C 1, i.e. diffusivity 1/f-noise is much weaker than its "seed" refernce level 
corresponding to D'/D = 1 . 

The opposite situations, when Tq/2Dto = D'/D ^ 1 , are possible too. Some examples can be 
found in [24j [25l HI] . Besides, recollect that in metals notably below Debye temperature, vice 
versa, velocity (direction) relaxation is slower than energy relaxation. 

This CF represents peculiar case of stable distributions, or the "a -stable Levy distributions", 
which corresponds to the "stability parameter" value a = 2 — , thus suggesting alternative to 
the Gaussian distribution for which a = 2 . 

Interestingly, this "quasi-Gaussian" [25] [55] distribution for the first time had appeared just in 
the commented preprint and independently, in same 1982, in mathematical work [60] (see also 
remarks in [28] ). In the standard classification of stable distributions (see e.g. [12]) such one 
is absent. In our considerations, it plays role of intermediate idealization between absolutely 
nonphysical Gaussian distribution and physically meaningful distribution characterized by (|2.33p 
or, more generally, (|2~2"4"]) - (j2~2"gj) . 

Qualitative difference between random walk processes represented by CF (|2 . 23[) and CF (|2.33[) 
(or (|2.24j) - (|2.28[) ) is that increments of the first of them are purely statistically independent and 
have infinite variance, while increments of the second have finite variance and all are essentially 
statistically dependent, regardless of their time separation. Thus, our cut-off of scale invariance 
( "fractality" [25]) of random walk r{t) from below (at small scales) completely changes its sta- 
tistical contents, especially in physical applications. 

Unfortunately, such treatment of "fractal random walks" is out of traditions of "mathematical 
physics" which develops in the firm belief that dynmical random walks in statistical mechanics 
always can be decomposed into "independent" fragments and thus reduced to stochastic random 
walks like e.g. "Levy flights" (perhaps, just by this reason the "Journal of Statistical Physics" 
resembles sooner a "book of problems" from probability theory than researches statistical me- 
chanics). However, as was predicted by Krylov [2B], that are vain hopes, first of all in respect 
to "molecular random walks" in many-particle systems [3T] [32] [M] . Applying the A.Einstein's 
sentence, "God does not play dice " ! 

R.D.Hughes, E.W.Montroll and M.F.Shlesinger, J.Stat.Phys. 28 111 (1982). 

This is key statement and key result of the preprint under consideration. Importantly, it appeals 
to rigorous statistical mechanics of many-particle systems, instead of probability-theoretical spec- 
ulations. Now, we can state that statistical mechanics indeed justifies this appeal (see [27] - [39], 
[4U1 l4"2l [H] and references therein) . 

As was shown in [25], probability weight of non-Gaussian (more or less long) tails of the r(t) 's 
distribution typically is around 0.03 . Nevertheless, this means that probabilities of large devi- 
ations of r(t) from typical diffusion rate (i.e. probabilities of \r(t)\ \j2Dt ) always are very 
mich greater than in Gaussian model (see [31] [37] [39] 148]). 
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Of course, statistics of large deviations is sensitive to details of the function F(z) , especially if 
it decays by a power law, e.g. F(z) = 1/(1 + z//3) 1+/3 (we took into account conditions (I2.27P - ). 
At that, however, to satisfy actual finiteness of all r(t) 's cumulants, we should multiply such 
function of z = r 2 /AD't b a cut-off factor F' {r 2 / v^t 2 ) , where F'(z') is fast decreasing function 
of z'r 2 jv^t 2 and vo is characteristic velocity of carriers (e.g. thermal or Fermi velocity). Then 
probability distribution of r(i) asymptotically (at v^t 2 / D't — > oo ) obeys scale invariance r 2 cx t 
but far tails of this distribution and its higher-order cimilant significantly deviate from r 2 cx t 
(therefore, in particular, the "noise of 1/f-noise" [25] can be of 1/f -type with 7 > 1 ). A similar 
picture arises in statistical-mechanical theory of "molecular random walks" [3pJ[3l] 32, 34, 39] , 148 ) . 

[63] Of course, that are not absolutely convincing reasonings, and one can suggest contrary instances 
to them. For example, in metals, due to the Fermi statistics of carriers, even under simple inelastic 
scattering (above the Debye temperature) r\ ss Ag S> 2Dtq IDt^ , so that the Hooge constant 
rises nearly Tp/T times as compared to the "seed" level [24] (as if only ~ (T/Tp)N of total 
number N of conduction electrons, - from vicinity of Fermi surface, - effectively contribute to 
1/f-noise). 

On some of other possible variants of relationships between characteristics micro-scales see [24] [25j 
[28] [301 EH @0j [41] . Unfortunately, exact equations of microscopic Hamiltoniam model of electron's 
random walk in phonon field [42] still are not analysed up to quantitative estimates to compare 
with the commented phenomenology. 

[64] Though, strictly speaking, this integrability is not consequence of stationarity in its physical 
sense. In other words, 1/f-noise can look formally non-stationary even in thermodynamically 
equilibrium (all the more, in steady non-equilibrium) systems! This important theoretical fact 
was demonstrated and explained already in [27, (se also [28 ] [30 ] ). 

[65] On this subject see also [271128] , 

[66] This seeming paradox is main (practically fatal!) obstacle in ways to adequate theoretical per- 
ception of 1/f-noise. 

Though, it is trivial logical truth that (long) statistical correlations can take place without (long) 
physical (cause-and-consequence) correlations beyond them. Moreover, forgetting of the past 
(history of a transport process) means impossibility to keep the future (fluctuations in transport 
rate) predictable and controllable. Therefore short-memory transport process inevitably acquires 
scaleless long-term fluctuations (and thus statistical correlations) of its instant rate. 

Nevertheless, physicists (even theoreticians) accept for a fact the prejudice that physical inde- 
pendence is the same as statistical independence. Consequently, they continue, like thirty years 
ago, to search sources of 1/f-noise observed in transport processes (and other irreversible and 
transition processes) everywhere except randomness of these processes in itself! 

Notice that measurements of effective "Hooge constant" ("Hooge parameter") in very different 
systems, - including modern small-size diamond, grathene and carbon nanotube based devices 
(see, for example, [57], [BS] and [5^1 ED], respectively, and references in [H]), - usually show values 
from rather narrow interval 10 -4 < a < 10 _1 , in spite of giant variety of conditions and regimes 
of charge transport under investigations. But, strikingly, this fact does not enforces investigators 
to assume that mysterious source of 1/f- noise is neither charge traps nor structural defects or 
impurities (or even phonon energy fluctuations, - see discussions in |29[ 142]). That is what the 
Koz'ma Pritkov's joke is about (see epigraph above). 

In reality, of course, the enumerated factors, along with many others, - e.g. Coulomb interaction 
between carriers and related screening and space charge effects, - all can essentially influence 
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1/f-noise, but this does not mean that they are its sources (similarly, they can influence white 
noise although not being its sources). 

Hence, the actual problem waiting for theoreticians is development of dynamical (i.e. statistical- 
mechanical) models of systems under practical interest (thus continuing our first steps made in 



The phenomenology we are commenting now is too primitive from viewpoint of such complicated 
inhomogeneous structures as mentioned various field-effect transistors [63 |68l ESI EH] • Neverthe- 
less, it can be notably helpful in interpretation of experimental data (see comment |78j below). 

[67] M.V. Hauf, L.H. Hess, J.Howgate, M. Dankerl, M. Stutzmann, and J. A. Garridoa, Appl. Phys. 
Lett. 97 093504 (2010). 

[68] M. Dankerl, M.V. Hauf, A. Lippert, L.H. Hess, S. Birner, I.D. Sharp, M. Stutzmann, J.A. Garrido, 

Adv. Funct. Mater. 20 3117 (2010). 
[69] J.Malnnik, I. Heller, A.M. Janssens, S.G.Lemay, and C.Dekker, Nano Lett. 8 (2) 685 (2008). 
[70] Ju Hee Back, Sunkook Kim, Saeed Mohammadi, and Moonsub Shim, Nano Lett. 8 (4) 1090 

(2008). 

[71] "Rate of diffusion" by its sense is a time-nonlocal quantity smeared over time intervals much 
longer than system's memory for diffusion. That is why this quantity has no definite instant 
values, is free of system's back reaction, and fluctuates with no certain time scale. 

[72] On this subject see also ^\ and discussions in [2"81 14"2"] . 

[73] In English, see also papers [15] plus [15] - [15] and other our works on the nonlinear FDR, or 
"generalized FDR" , referenced in [46] - [48] . 

[74] See additionally [24] and references therein. Some usefull information can be found also in [54) . 

[75] Though, see above comment (64! . Here, we again have to distinguish between formal and physical 
stationarities, since they are not equivalent, and integral in (|3.12[) can be unboundedly growing 
function of t although being characteristics of equilibrium state. This means [27] that result of 
conductance measurement depends on its duration t but in no way on time when it begins, and 
at that does not improve when t increases (in opposite, longer time averaging makes it worse, 
since increases its mean-square relative deviation from ensemble average!). In such case (|3.12p 
must be treated like the structure function of non-stationary random process (for examples see 
|27[l30j). and relations between theoretical expressions and experimental procedures need in more 
careful analysis. 

[76] To be precise, it anyway decreases with time faster than <x 1/t, while K g (t) by its sense decreases 
slower than oc 1/t (if not increases, - see above). 

[77] Clearly, the CF S(iu\x) here is the same as A 00 (m|x) in absence of low- frequency flicker fluc- 
tuations which are separated into factor £(t) . At that, the latter can be treated as derivative of 
randomly non-uniform "inner time" of the transport process, J £(i') dt' , in respect to real time. 

Similar separation of high-frequency noise and 1/f-type noise appears also, in the course of nat- 
uralal approximations, in the microscopic statistical-mechanical theory of molecular Brownian 
motion (see e.g. [301 ED IM1 IB] ) ■ 

At the same time, importantly, formula ()3.20j) can be derived without such separation, by using 
the equilibrium scale invariance reasonings plus some FDR only (see |24j). 



[Zl], [27] - [42], @8]). 
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[78] Though, after introducing effective number of carriers under simultaneous transitions through 
the junction, N(x) = \ J{x)\ro/e , and rewriting e\J(x)\/ro = J (x)/N(x), formula (|3.22p looks 
like in Ohmic regime. 

The general formula (|3.20[) covers cases when, firstly, the conductor is neither discrete junction 
nor homogeneous distributed sample and, secondly, charge carriers can interact one with another, 
that is be (physically and statistically) mutually dependent. 

In particular, we may apply formula (|3.20[) to field-effect transistors (FET), as far as it is possible 
to neglect their gate currents and treat FET like two-terminal conductor with variable properties 
tuned by gate voltage V g . At that, according to experimental observations (see e.g. [57J - [70]). 

2 _______ 

usually Sj oc J at fixed V g . Therefore, firstly, p. 201) can be written as 

At \ = 111 - 1a( x ) a o 
1 ' ~ J 2 S(x)r (x) ' 

where now x — V g — V t h is operating gate voltage counted from some "threshold voltage" Vth , 
f = u/2ir, ao = (l/3)(ln |u;|to)~ 2 ~ 0.001 is the seed Hooge parameter, and S(x) as before 
is white noise power spectral density (of source-drain current). Secondly, we can use the Nyquist 
formula, S(x) = 2TG{x) , with G(x) being source-drain conductance. 

For instance, consider charge transport in single-wall carbon nanotube (SWCNT) based FET. If 
length L of SWCNT and x are large enough, then G G = (4e 2 /2ttH) X/L [75], where A , as 
above in Sec. 2. 7, is mean free path. As there, we have to take go = er a/L and naturally assume 
that r — X and r = X/vp , where Vp is Fermi velocity. Then 

ttHvf a e 2 

A = ~4TL~ ' 00 = N = 2TC~ q ' 00 ' 

where N = (2/ttHvp) e\x\L is number of carriers (with 2/ttHvf representing ID density of 
states in "hyperbolic" or"conical" band), a= (e\x\/2T) ao (with e\x\ representing Fermi energy, 
so that ratio a/a looks like for metals), and C q — e 2 (2/ttHvf) L is SWCNT's "quantum 
capacitance" . 

For numerical estimates, we must take into account that vf ~ 10 8 cm/s [79], and, - in view of 
just made assumptions, - consider comparatively long tubes with L ~ 3 /im [69, 70 at sufficiently 
large \x\ > 0.3 V (see e.g. Fig.l in [7D] for visuality). Then the above formula yields A ~ 40 -5 . 
It agrees with A measured in [55J, as well as with minimum (measured under \x\ > 0.3 V) of A 
values found in 70 . At that, V ~ 2000 and a w 10oo ~ 0.01 , in good agrreement with values 
a ~ 0.01 typically observed in [70] . 

At |x| < 0.3 V, in [70] significantly greater A and a values were obtained. Apparently, they have 
subjective reason: the number of carriers there was estimated like we have done above, by N — 
C q \x\/e = No(x) , although in reality at relatively small |x| (when metallic 'hyperbolic-band" 
conductance through SWCNT channel still is not formed) N = N(x) decreases in a more fast way. 
Namely, such that N(x)/N (x) ~ G{x)/Gq <C 4 at small |x| , as simple analytical consideration 
does show (at that G(x) cx J changes like current at Fig.l in [70], and A^ 1 (x) oc G(x)). By 
this reason, a appears to be strongly overestimated. Besides, objectively, high a values may be 
attributed to inhpmogeneity of gating the channel at small \x\ (since any inhpmogeneity lowers 
conductance and rises effective Hooge parameter). 

Analogous satisfactory and even better situations, - as to agreement of our theory with experi- 
ment, - take place in cases of FET exploiting "H-terminated diamond" and graphene micrometer- 
size sheets [571155]. All that will be considered elsewhere. 
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[79] M.J. Biercuk, Sh.Ilani, Ch.M. Marcus, and P.L.McEuen, in [A.Jorio, G.Dresselhaus, 
M.S.Dresselhaus (Eds.). Carbon Nanotubes. Topics Appl. Physics 111, 455493 (2008), Springer- 
Verlag, 2008]. 

[80] It is tercentennial heritage of the "law of large numbers" mathematically deduced by J. Bernoully 
[8Tj from assumption that seemingly (physically) independent observations (events) are statisti- 
cally independent (in the sense of modern probability theory). This assumption allowed Bernoulli 
to manage with a small set of particular a priori "probabilities of events" and to propose their 
determination by means of time averaging. Indded, otherwise it would be necessary to introduce 
probabilities of pairs, ternaries, etc. (up to infinity), of events, so that number of ndeterminates 
always would stay greater than number of observations. 

However, in the twentieth century N. Krylov 26 showed that in statistical mechanics such as- 
sumptin in general is wrong, that is "probabilities of events" (relative frequencies of events) 
significantly vary from one phase trajectory (experiment) to another, regardless of duration of 
timr averaging. Thus, the "law of large numbers" does not work, and statistics of real "flows of 
events" is essentially non-Gaussian even in low-frequency domain (see also remarks in introduc- 
tion sections of [30 | 134 ) [35 ] ) . 

Therefore, A. Kolmogorov [52] was right recommending scientists be careful with assumptions 
about "independency" of physical random events. 

At this point, it is time to follow the Koz'ma Prutkov's advice: "If you have a fountain, stop it 
up. Let it too have a rest". 

[81] Jacobi Bernoulli. Ars conjectandi. Basel, 1713 [in English: Jacob Bernoulli and E.D.Sylla (trans- 
lator). Art of conjecturing. John Hopkins Univ., 2005; in Russian: Ya. Bernulli. O zakone bol'shih 
chisel. Moscow, Nauka, 1986]. 

[82] A.N. Kolmogorov. Foundations of the theory of probability. Chelsea, New York, 1956 (Russian 
original first was published in 1936). 
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